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Introduction 


This  report  summarizes  work  carried  out  under  grant  number  AF- AP0SR-719"65 • 
Most  of  the  effort  was  devoted  to  obtaining  a  better  understanding  of  the  be¬ 
havior  of  cuprous  oxide  at  high  temperatures  with  special  reference  to  the 
electrical  properties  and  defect  structure.  In  this  work  we  were  lead  to  an 
investigation  of  other  transition  metal  oxides.  This  phase  of  the  work  is  being 
continued  under  grant  number  AFOSR-68-1371  so  that  in  some  respects  this  docu- 
ment  will  be  in  the  nature  of  a  progress  report.  This  is  true  also'because  a 

■  V  .) 

considerable  amount  of  the  early  effort  was  devoted  to  construction  of  apparatus 
that  will  be  used  in  subsequent  work  under  the  new  grant. 

In  subsequent  sections  we  summarize  the  main  achievements  of  the  work. 
Included  as  an  appendix  are 

•*  (a)  Diffusion  in  Oxides  and  Sulfides  pi  rented  at  the  International 

Symposium  on  Sintering  and  Related  Phenomena  at  Notre  lame  university.  As  the 
publication  of  the  proceedings  of  this  symposium  has  been  delayed  the  manu¬ 
script  is  included  nere  for  wider  circulation. 

(b)  Preprints  of  two  papers  Just  completed  and  to  be  published  as  part  of 
the  Proceedings  of  the  Symposium  on  Transport  in  Oxides  held  at  N.B.S.  this 
month  (October  1967). 

Studies  on  Cuprous  Oxide 

The  major  effort  here  involved  construction  of  an  A.C.  Ball  effect  apparatus 
for  measurements  on  low  mobility  materials  at  high  temperatures.  The  main  re¬ 
sults  of  measurements  on  pure  crystals  have  been  published  -  ".High  Temperature 
Hall  Effect  in  Cuprous  Oxide”  by  H.  L.  McKinzie  and  M.  O'Keeffe  and  will  not 
be  discussed  farther  here. 

Two  types  of  thermogravimetric  measurement  have  been  made  with  cuprous 
oxide.  In  the  first  the  kinetics  of  vaporization  were  studied  in  various 


atmospheres.  TLe  results  of  this  work  were  presented  at  the  International 
Symposium  on  Kindles  of  Ionic  Reactions  at  Alfred,  N.Y.  and  will  appear  in 
the  published  proceedings. 

Attempts  were  elso  made  to  study  the  kinetics  of  change  of  stoichiometry 
of  the  oxide  following  a  rapid  change  of  composition  of  the  surrounding  at¬ 
mosphere.  The  rapid  rate  of  weight  loss  by  evaporation  precluded  meaningful 
measurements  so  *n  electrical  method  was  developed  for  making  these  studies. 

The  results  of  this  work  are  included  in  this  report  as  an  appendix. 

Studies  with  other  Oxides 

The  lnterdiffuslon  studies  by  thermogravimetric  means  have  been  extended 
to  wdstite  (FeOx).  Here  there  is  considerable  interest  in  the  effect  of  defect 
interactions  in  this  very  non- stoichiometric  compound.  Again  a  preliminary 
account  of  this  work  is  included  with  this  report. 

Finally  an  extensive  set  of  measurements  have  been  made  on  the  thermo¬ 
electric  power  and  conductivity  of  manganous  oxide  (M 10).  To  a  certain  extent 
these  measurements  have  been  overlapped  by  the  work  of  Hed  and  Tannhauser 
[J.  Chan.  Phye.  kj,  2090  (1967)].  However  there  are  still  some  puzzling  problenu 
in  the  interpretation  of  these  results  and  their  analysis  is  not  yet  complete. 

It  is  hoped  that  experiments  now  underway  with  doped  crystals  will  shed  more 
light  on  this  problem. 
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Over  the  last  ten  or  fifteen  years  a  large  tody  of  experimental  data 
has  been  accumulating  on  diffusion  in  oxides  and,  to  a  lesser  extent,  in 
sulfides*  In  this  review  these  data  are  collected  together,  and  an  attempt 
is  made  to  indicate  the  more  interesting  of  recent  developments. 

In  contrast  to,  say,  the  alkali  halides,  the  oxides  present  an  enormous 
variety  of  materials  with  very  different  properties  so  that  it  is  difficult 
to  treat  the  subject  in  general  terms.  Because  of  this  diversity  of  behavior 
too,  there  has  been  little  tendency  to  phudy  any  one  system  in  depth,  so 
that  despite  the  number  of  systems  that  have  been  studied,  it  is  probably 
true  that  in  no  instance  is  our  knowledgeAthe  nature  and  mobility  of  de¬ 
fects  as  complete  as  it  is  for  example  in  KC1  or  AgCl. 

In  an  attempt  to  systematise  the  discussion,  it  is  convenient  at  the 
outset  to  divide  oxides  and  sulfides  into  three  broad  types  as  follows: 

A.  Materials  with  a  very  small  concentration  of  native  defects.  v  Very 
small'  may  be  taken  to  mean  less  than  that  detectable  by  conventional 

'e 

t 

chemical  means  (say  <  .01$).  Uhder  this  classification  will  come  the 
oxides  and  (to  a  lesser  extent)  the  sulfides  of  the  metals  of  the  first 
three  groups  of  the  periodic  table.  LigO,  MgO,  AlgOy  MgS,  etc. 

B.  Materials  in  which  the  defect  concentration  is  more  readily  measurable, 
but  in  vhich  the  defects  may  be  considered  to  be  in  dilute  solution 
(say  .01  to  1$),  These  materials  will  be  largely  transition  metal 
compounds  In  vhich  more  than  one  cation  valency  is  possible.  Ex¬ 
amples  are:  CoO,  N10,  Cu^O. 

C.  Materials  vith  a  large  (>  1$)  concentration  of  defects.  Well  known 
examples  are  TiO,  FeO,  FeS. 

The  distinction  between  types  B  and  C  has  recently  been  elaborated  on 
by  Anderson. ^  It  is  to  be  understood  that  the  classification  above  is 
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very  rough  and  ready  and  that  a  given  material  may  veil  fall  into  different 
types  according  to  conditions  of  temperature  and  ambient  atmosphere. 

We  will  now  discuss  diffusion  prob3.ems  that  are  pev.  liar  to  each  type 
separately. 

Type  A 

Materials  falling  under  this  heading  will  have  a  very  low  concentration 
of  native  defects,  which  may  therefore  he  expected  to  show  nearly  ideal  be¬ 
havior,  and  there  is  little  doubt  that  diffusion  in  a  monocrystal  will  occur 
by  a  defect  mechanism. 

However,  for  this  type  of  oxide  there  is  a  large  experimental  obstacle 
to  the  understanding  of  the  defect  structure  and  hence  of  diffusion  mech¬ 
anisms.  This  is,  of  course,  the  problem  of  obtaining  crystals  of  sufficient 
purity  that  intrinsic  properties  of  the  crystal  can  be  studied.  For  mate¬ 
rials  of  this  type  the  energy  of  formation  of  a  defect  is  very  large.  If 

ve  accept  2ev  es  a  plausible  value  for  the  energy  of  formation  of  a  de- 

(2) 

feet,'  '  then  the  impurity  concentration,  c,  expressed  as  a  fraction  of 
total  sites  will  be 

-8 

c  ~  exp  (-2ev/kT)  ~  10  at  1000 °C 

It  will  require  only  e  very  small  concentration  of  impurities  (e.g.  Fe^ 
in  MgO  or  Al  in  TiOg)  to  produce  defect  concentrations  greatly  in  excess 
of  this  intrinsic  concentration  in  the  pure  crystal.  It  is  very  likely 
then,  that  defect  controlled  properties  such  as  self  diffusion  will  be 
determined  by  the  fortuitous  defect  concentration.  This  is  made  clear  by 
the  work  of  Wert*  and  his  collaborators  who  have  made  er.  ext^^sive  etudy 
of  defects  in  MgO  and  concluded  in  'To  our  knovlvlge,  no  MgO 

crystals  have  ever  been  grown  which  did  not  show  ar.  ESR  absorption  by  one 
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or  more  impurity  ions.  In  the  best  specimen  we  have  observed  thus  far,  iron 
was  present  to  the  extent  of  5  parts  per  million.  It  is  likely,  therefore, 
that  many  of  the  electronic  properties  of  magnesium  oxide  measured  to  date 
are  to  a  serious  degree  affected  by  the  presence  of  impurities.  In  addi¬ 
tion  to  direct  effects  of  the  ions  themselves,  one  has  indirect  effects 
arising  from  positive-ion  vacancies  ascribable  to  the  presence  of  trivalent 
ions.' 

In  the  some  vein,  attention  is  drawn  to  the  work  of  Kanzaki  and  co- 
(k\ 

workers'  '  who  notice  significant  effects  of  impurities  on  ionic  conductivity 
of  KC1  with  divalent  impurity  concentrations  as  low  as  0.01  p.p.m. 

Despite  these  rather  pessimistic  remarks,  it  is  worth  enquiring  in  gen¬ 
eral  terms  what  type  of  behavior  is  expected  of  oxides  of  different  struc¬ 
tural  types.  Let  us  consider  for  example  some  simple  oxides  of  general 
formula  M^O,  MO  and  MO^  and  enquire  as  to  whether  anion  or  cation  mobility 
will  be  the  greater. 

If  the  crystals  are  considered  to  be  an  array  of  charged  ions,  it  is 
fairly  clear  that  a  very  important  factor  in  determining  the  type  and  mo¬ 
bility  of  defect  in  the  crystal  will  be  the  variation  of  the  electrostatic 
potential  (the  Madelung  potential)  from  point  to  point  in  the  crystal. 

The  electrostatic  potential  at  an  ion  site  in  a  simple  ionic  crystal 
may  readily  oe  evaluated  in  terns  of  certain  lattice  sums^.  Values  that 
are  typical  are  given  for  three  oxides  in  Table  1.  In  the  table  $c  and 
are  the  electrostatic  potentials  at  a  cation  and  anion  site  respectively 
arising  from  the  presence  cf  all  the  other  ions  in  the  crystal. 

Die  significance  of  these  potentials  is  as  follows.  The  work  necessary 
to  remove  an  ion  from  the  otherwise  perfect  crystal  is  -2#  where  Z  is  the 
charge  on  the  ion.  2$  varies  roughly  as  2  (cf  Table  1).  As  a  result  ve 


Table  1 


Electrostatic  potentials  at  ion  sites  in  oxides 


Oxide 

Structure 

$c  (volts) 

*A  (volts) 

-y>v) 

•Va(cV) 

Na20 

anti-fluorite 

-10.6 

•  19. -2 

10.6 

38.4 

MgO 

rock  salt 

-25.9 

23.9 

1*7.8 

1*7.8 

Zr02 

fluorite 

-1*3.1 

23.2 

172.1* 

1*6.1* 
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might  expect  that  the  formation  of  defects  in  which  cations  are  absent  from 
their  regular  sites,  will  most  readily  occur  when  the  cation  charge  is  small 
Alternatively,  if  cation  and  anion  vacancies  are  present  in  equal  numbers 
(Schottky  defects  in  a  stoichiometric  crystal)  one  would  expect  cation  mo¬ 
bility  to  be  greatest  when  the  cation  charge  is  small. 

In  the  actval  energetics  of  the  formation  and  migration  of  defects  in 
ionic  crystals,  relaxation  (i.e.  polarization)  of  the  crystal  around  the 
defect  will  be  very  important.  In  general  polarization  will  favor  the  for¬ 
mation  and  migration  of  vacancies  of  the  less  polarizable  ion^  (normally 
this  is  the  cation)  so  that  in  MO  (and  possibly  even  in  as  well  as 

rcr  MgO,  cation  diffusion  is  likely  to  be  more  rapid  than  anion  diffusion. 
Conversely  for  MOg  and  higher  oxides  it  is  expected  that  anion  diffusion 
will  be  more  rapid. 

Extension  of  this  type  of  argument  suggests  that  the  same  general  rule 
will  also  apply  to  non-, stoichiometric  crystals. 

It  should  be  remarked  that,  given  the  problem  of  obtaining  sufficiently 
pur?  crystals,  a  profitable  course  to  pursue  with  this  type  of  oxide  would 
be  to  deliberately  dope  the  crystal  with  altervalent  impurities.  Ifcere 
appears  to  have  been  very  little  work  done  along  these  lines  (with  the  ex¬ 
ception  of  work,  on  certain  fluorite- type  oxides  mentioned  below)  in  con¬ 
trast  to  the  very  extensive  work  reported  on  doped  alkali  halides.  It  would 
be  very  interesting  to  study,  for  example,  MgO  doped  with  small  quantities 
of  AlgOj. 

For  this  type  of  oxide  (type  A)  there  is  the  attractive  feature  that 
electronic  conductivity  la  sometimes  negligibly  small  compared  with  ionic 
conductivity.  It  is  usually  true  also  (as  discussed  above)  that  one  ion 
la  considerably  more  mobile  than  the  other.  If  this  is  the  case,  the 
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dif fusion  coefficient  D,  for  the  most  mobile  species  is  readily  obtained 
from  the  conductivity  o  using  the  Nemst-Einstein  equation^ 

D/a  -  pcT/CNq^  (1) 


Here  N  is  the  number  per  unit  volume  of  ions  of  charge  q.  ^  is  a  n>imerical 
factor  -  the  correlation  factor  -  introduced  to  take  into  account  that  suc¬ 
cessive  Julies  of  the  defect  may  not  be  entirely  random.  Normally  J  is  close 
to  unity. ^ 

It  must,  of  course,  be  established  that  the  conductivity  is  ionic.  This 
has  been  done  for  remarkably  few  oxides.  Transport  numbers  can,  however,  be 
measured  fairly  readily  from  e.m.f.  measurements.  Consider  for  example  a 
concentration  cell  of  the  typ2 


Pt.  0^) 


oxide 


Pt,  o2(p£) 


If  the  oxide  were  an  ionic  conductor  the  e.m.f.  of  the  cell  would  be 
given  by  the  familiar  Nerast  equation: 


■Wc  ’  *  <VP2)V2 


(2) 


If,  on  the  other  hand,  there  were  an  appreciably  electronic  conduction, 
this  would  have  the  effect  of  partially  'short  circuiting'  the  cell  and  the 

/g\ 

e.m.f.  would  be  reduced  to  '  ’ 


K 


ti  Sionie 


(» 


(a ) 

where  is  the  ionic  transport  number.  Sciraalsried'"  bss  recently  given 

s  detailed  discussion  of  this  type  of  measurement  and  results  for  several 

% 

oxides. 

As  an  example  of  the  utility  of  electrical  measurements  of  this  type 
ve  say  take  as  an  example  the  results  obtained  co  the  mixed  oxide  Zr^^Ca^O,,^. 
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For  this  oxide  the  ionic  transport  number  is  unity,  and  the  oxygen  ions  are 
the  most  mobile  species.  D  for  oxygen  has  been  measured  directly  by  tracer 
methods  (see  Table  3)  and  is  in  excellent  agreement  with  that  calculated 
from  electrical  conductivity  using  equation  (l),^^  The  importance  of  trans¬ 
port  number  measurements  is  illustrated  by  results  for  ThCL-YnCL  solid  so- 

2  c.  3 

lution3  which  might  be  expected  to  behave  very  similarly  to  ZrOg-CaO,  but 
which  in  fact  showed  predominantly  electronic  conductivity.^^  Similar 

fl2) 

measurements  have  also  been  made  on  pure  MgO  but  with  conflicting  recults.  ' 
Type  B 

For  this  type  of  oxide  the  situation  is  rather  different.  Here  devi¬ 
ations  from  stoichiometry  are  sufficiently  large  that  the  native  defect  con¬ 
centration  is  far  greater  than  the  normal  level  of  foreign  atom  concentration 
in  the  absence  of  deliberate  doping.  In  these  materials  too,  the  conductiv¬ 
ity  is  almost  certainly  electronic. 

lii  several  instances  it  has  proved  possible  to  measure  departures  from 
stoichiometry,  either  by  chemical  means  (e.g.  by  analysis  for  Co^  in 
Co0^+x^')  or  gravimetrically  (e.g,  for  CUgO^  If  it  is  assumed 

that  one  type  of  ionic  defect  predominates  in  the  material,  the  analysis 
will  give  the  defect  concentration  directly.  As  an  example,  if  the  only 
ionic  defect  in  Co0^+^  were  cation  vacancies  (as  is  thought  to  be  the  case), 
then  there  are  x  vacancies  per  cobalt  ion. 

Ih  general  the  diffusion  coefficient  D  is  given  in  terms  of  the  atom 
fraction  c^  of  defects  responsible  for  diffusion  by 

D  -/  Vi 

Here  is  the  diffusion  coefficient  for  the  defect  and  f  is  the  correlation 
factor  mentioned  already.  Knowledge  of  c^  therefore  enables  to  be 
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calculatei  from  the  tracer  diffusion  coefficient  D. 

For  a  dilute  solution  of  defects  we  can  express  the  temperature  and 
pressure  dependence  of  c^  in  terms  of  thermodynamic  quantities  as  follows. 

We  take  as  an  example  an  oxide  MO  in  which  cation  vacancies  and  M  ions  are 
formed  in  incorporating  an  excess  of  oxygen.  In  terms  of  Kroger’s  notation 
the  reaction  is 

1/2  02  — »  0Q  +  Vj»  +  h#  (5) 


Here  0Q  is  an  oxygen  ion  on  an  oxygen  ion  site.  is  a  vacancy  made  by  re¬ 
moving  a  M  atom  from  a  metal  ion  site.  We  have  assumed  that  this  vacancy  is 
’ionized’,  i.e.  that  it  has  lost  one  of  the  positive  charges  associated  with 
this  site.  The  positive  charge  is  represented  by  h®  -  it  may  alternatively 
be  thought  of  as  a  M  ion  in  the  crystal.  The  equilibrium  constant  for  the 
reaction  (5)  is 


K  *»  exp 


*> 


a(h°)  a(v^)  a(0Q) 
[  a(02)]  ^ 


(6) 


The  activity  of  0Q  is  constant  and  may  be  chosen  to  be  unity.  For  a 
dilute  solution,  we  may  replace  the  activity  of  defects  by  concentrations: 

a(V^)  «  a(h®)  »  (7) 

and 

»(02)  -  P0  (8) 

*  °2 

Prom  (6),  (7)  and  (8)  ve  have  immediately 

c  *-  exp  (-4G®/2RT)  P  ^  (9) 

Note  that  LG*  is  the  standard  free  energy  change  for  reaction  (5). 

At  constant  oxygen  pressure  ve  have 
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5  In  «  -aH4 

FPTtJ  ~ 


do) 


AH*  Is  now  the  standard  enthalpy  change  of  reaction  (5),  or  stated  differ¬ 
ently  AH*  is  the  heat  of  solution  of  1  g.  atom  of  oxygen  in  M0^+J£. 

Note  that _ if >  instead  of  reaction  (5)  ve  had  the  vacancy  doubly  ionized: 


1/2  02  — >  0o  +  Vj>  +  2h* 


(11) 


then  would  be  given  by 


Cj  f  (1/4)1/5  exp  (-AQ°/5RT)Po  (12) 


In  general 

e,  &  exp  (-AGVnRT1)  ?  1/2n  (13) 

1  o2 

and 


atoci  -  ^ 

nK 


(14) 


where,  as  before.  Ah*  is  the  heat  of  solution  of  1  g.  atom  of  oxygen  in  the 
oxide.  The  value  of  n  can  be  determined  from  the  oxygen  pressure  dependence 
cf  cA  (equation  (13)). 

Dj  will  likewise  depend  or.  temperature.  A  convenient  expression  for 
is  given  by  the  transition  state  theory  of  Eyring.^^ 

Da  -  (kT/h)  a2  exp(AsVR)  exp  (-aH*/?T)  (15) 

AS*  and  AH^  may  be  thought  of  as  respectively  the  entropy  and  enthalpy 
cf  activation  for  the  elementary  Jump  process  in  diffusion,  a  is  the  Jump 
distance  and  k  and  h  are  Boltzmann's  and  Planck's  constants  respectively. 

From  equation  (4),  (13)  and  (15)  It  is  clear  that  if  and  D  are  known 
as  functions  of  temperature,  AH*,  AH*,  AS*  and  <\S*  may  be  found. 
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Data  for  several  oxides  have  been  analyzed  by  O'Keeffe  and  Moore'  . 
Table  2  lists  AH$>  AH0,  AS^  and  AS0  for  oxides  in  which  cation  vacancies 
are  supposed  to  be  the  major  defect.  Note  that  in  the  table  AH°  and  AS0 
refer  to  the  formation  of  1  mole  of  vacancies  so  that  in  the  case  of  CUgO 
AH0  is  one  half  of  the  heat  of  solution  of  a  gram  atom  of  oxygen.  Shis  . 
arises  quite  simply  as  in  Cu^O  there  ere  two  cation  sites  for  each  oxygen 
site. 

Referring  to  Table  2,  it  can  be  seen  that  AS”^  is  small.  This  is  to  be 
expected  for  diffusion  by  a  point  defect  mechanism  involving  very  few  atoms 
in  the  elementary  jump  process.  For  the  isomorphous  series  of  monoxides, 

AH^  varies  only  very  little;  this  is  a  valuable  confirmation  of  the  correct¬ 
ness  of  the  assumption  that  diffusion  is  by  the  same  mechanism  in  all  these 
oxides. 

Turning  now  to  the  quantities  for  formation  of  a  vacancy,  it  may  be 
* 

seen  that,  in  contrast  to  AHt,  aH°  varies  very  dramatically,  and  it  is 
instructive  to  see  why  this  should  be  so. 

We  consider  vacancies  to  be  formed  by  the  overall  reaction  (5)  (this 
is  confirmed  by  the  oxygen  pressure  dependence  of  c^  found  for  CoO).  Now 
in  terms  of  the  'real1  charges  in  the  crystal  M^ 02~  reaction  (5)  is 

<  *  V2  °2  -•*  0*-  +  v;  *  # 


which  may  be  broken  down  into  the  following  steps: 


<•>  if- 

(b)  2e”  4*  1/2  Og 

(c)  M2*  4  o2 

ges  gas 


vu  +  M2* 

M  gas 


rt2' 

0  gas 


$ 


+  0 


2‘ 


(5*) 


Table  2 


Thermodynamic  quantities  for  formation  and  migration  of  defects  in  several  oxides 


Oxide 

AH’ 

AS* 

AH* 

AS* 

Mn0a 

0 

FeO 

50 

-28 

CoO 

27 

-0.8 

10 

-12 

NiO 

23 

(-D 

ItO 

(-14) 

CUgO 

12 

-1.5 

22 

4.9 

Note:  Enthalpies  in  kcal  mole’1,  Etatropies  in  cal  deg’1  mole’1. 
a  reference  (lOO) 
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<d>  VM->VM  +  •' 

M  $  *  e' 

For  the  closely  related  series  of  oxides  MnO,  CoO,  FeO,  NiO,  the  energy 
involved  in  each  of  the  steps  (a)  -  (d)  would  be  ejected  to  vary  only  slightly, 
(to  the  other  hand  step  (e)  is  the  ionization  of  and  the  energy  for  this 
step  will  vary  considerably  from  atom  to  atom  (ranging  from  30.6  eV  for  Fe 
to  36.2  eV  for  Ni).  It  is  expected,  therefore,  that  the  variation  in  dH° 
will  be  mirrored  in  the  variation  of  the  third  ionization  potential  of  the 
cation.  Figure  1  shows  that  this  is  indeed  the  case.  It  can  be  seen  that 
this  type  of  non-stoichiometry  is  unlikely  for  MgO,  as  the  third  ionization 

ft 

energy  of  Mg  is  very  high  (80.1  eV)  . 

# 

Because  of  this  high  ionization  energy,  one  would  expect  in  fact  that 
a  hole  in  MgO  would  be  localized  on  an  oxygen  ion  (rather  than  on  Ms  ). 

This  appears  to  be  observed  in  electron  spin  resonance. 


Finally  it  should  be  remarked  in  this  connection  that  the  evaluation  of 

requires  the  assumption  that  one  type  of  defect  predominates  in  the  ma- 
(17) 

terial.  However  Both'  '  has  shown  that  in  FeO  there  is  an  appreciable  con¬ 
centration  of  interstial  iron,  so  that  reaction  (5)  should  possibly  be  replaced 
by  a  reaction  such  as 

1/2  02  — >  0Q  +  Fe"  +  2Vje  +  2h*  (l6) 

How  far  this  observation  affects  the  interpretation  of  results  for 
oxides  such  as  CoO,  RiO  and  Cu^O  (in  all  of  which  the  predominant  defect  is 
assumed  to  be  cation  vacancies)  is  not  known.  It  is  worth  noting,  however. 
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that  for  molten  salts  of  the  alkali  halides  there  is  good  evidence  for  inter¬ 
stitial  (in  the  sense  of  tetrahedrally  coordinated)  cations. 

It  should  also  he  mentioned  that  information  about  may  also  he  ob¬ 
tained  directly,  from  the  kinetics  of  annealing  of  a  crystal  containing  a 

non- equilibrium  concentration  of  defects.  An  example  is  to  be  found  in  the 
(to) 

work  of  Thomas,  '  who  measured  the  decay  of  conductivity  of  ZnO  contain¬ 
ing  an  excess  of  interstitial  zinc.  Zinc  is  supposed  to  dissolve  in  zinc 
oxide  according  to 

Zn  — >  ZnJ  +  e*  (17) 

so  that  the  conductivity  is  proportional  to  the  concentration  of  excess  zinc 
in  the  crystal  and  the  rate  of  change  of  conductivity  is  proportional  to  the 
rate  of  diffusion  of  Zn*. 

This  type  of  experiment  involves  diffusion  in  a  concentration  gradient 
(i.e.  interdiffusion)  -  a  topic  to  vhich  we  return  below. 

'JXEZS 

A  large  number  of  oxides  and  sulfides  exhibit  a  wide  range  of  non- 
stoichiometry.  Formerly,  thi*  was  described  in  terms  of  point  defects  as 
for  materials  of  type  B.  Hcvever,  detailed  studies  of  individual  systems 
have  frequently  revealed  the  existence  or  r  number  of  closely  related 
phases  of  definite  composition.  Particularly  notable  are  the  x-ray  studies 
of  Magneli  and  his  school  and  the  phase  equilibria  studies  of  fyring  and 
co-workers.  Examples  of  such  phases  are:  TiO^^  with  n  -  4-10  known, 
and  and  Mo0^_x  which  form  a  large  number  of  phases  of  composition 

8X1(1  \°3n-l*^21^  extreme  example  is  Bare 

earth  oxides  in  the  range  MgO^  -  MOg  afforded  another  example  in  which  the 
existence  of  a  number  of  stable  phases  of  definite  stoichiometry  2 
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(ox) 

has  been  established.' 

These  results  are  not  surprising  as  appears  on  first  sight.  Some  time 
(°h) 

ago  Bertaut  pointed  out  that,  from  considerations  of  electrostatic  energy, 
a  structure  with  a  large  concentration  of  vacancies  would  be  expected  to  be 
very  unstable  with  respect  to  an  ordered  phase.  Indeed  the  possibility  of 
existence  of  a  very  non- stoichiometric  compound  has  recently  been  questioned. ^ 

It  is  profitable  to  examine  the  situation  in  a  little  more  detail.  Pig- 
tire  2a  shows  a  two-dimensional  crystal  of  nominal  composition  AB.  A  tew  A 
atoms  have  been  removed^  and  clearly  those  empty  A  sites  should  now  be 
referred  to  as  A  vacancies.  Suppose  now  that  50$  of  the  A  atoms  wore  re¬ 
moved  in  a  systematic  way  giving  a  compound  ABg.  This  process  has  been  al¬ 
most  completed  in  Figure  2b;  it  can  be  seen  that  some  of  the  sites  that  were 
A  sites  in  AB  are  interstitial  sites  in  AB£.  So  far  there  is  no  difficulty. 

Let  us.  now  make  a  solid  solution  of  AB  and  ABg  such  as  shown  in  Fig¬ 
ure  2c,  It  Is  not  clear  any  longer  whether  we  should  refer  to  this  as  AB 
with  A  vacancies  or  AB,,  with  interstitial  A.  However,  as  mentioned  above, 
it  Is  in  fact  rather  unlikely  that  we  should  have  a  random  solid  solution 
of  this  type;  the  reasoning  is  as  follows. 

Referring  again  to  Figure  2b,  it  can  be  seen  that  each  interstitial  A 
ion  has  four  nearest  neighbor  A  ions — a  configuration  that  might  be  expected 
to  be  rather  unfavorable  for  electrostatic  reasons.  Magneli^2^  has  pointed 
out  that  the  same  stoichiometry  can  be  maintained  and  yet  have  no  more  than 
two  nearest  neighbor  A  ions  for  any  given  A  ion.  This  can  be  done  as  shown 
in  Figure  3a  which  is  taken  from  Magneli's  paper.  The  feature  to  note 
new  Is  that  most  of  the  ions  %  in  the  crystal  have  the  same  immediate  environ¬ 
ment  as  in  pure  ABg;  the  exceptional  ions  are  arranged  along  tig- rag  lines 
that  run  through  the  crystal.  Thus,  there  are  regions  of  pure  AB0  separated 
by  randomly  aligned  regions  of  different  stoichiometry  (richer  in  A).  - 
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Adjacent  ABg  regions  are  no  longer  in  phase,  and  to  this  extent  there  is 
strictly  no  longer  long  range  order  in  the  crystal.  However,  if  the  perfect 
ABg  regions  are  sufficiently  large,  the  x-ray  diffraction  pattern  of  the  non- 
stoichiametric  compound  will  be  that  of  ABg. 

We  enquire  now  as  to  ionic  mobility  in  a  crystal  of  this  type.  It  is 
assumed  that  it  is  a  very  unfavorable  configuration  when  an  A  has  more  than 
two  A  nearest  neighbors,  so  that  Jumps  that  result  in  this  type  of  configura¬ 
tion  will  not  occur.  This  rule  severely  limits  the  number  of  A  ions  that  can 
move  at  any  one  time.  Figure  5b  has  been  obtained  by  moving  two  ions  that 
can  move  with  the  above  restriction;  it  can  be  seen  that  in  the  new  arrange¬ 
ment,  ions  that  were  previously  unable  to  move  can  now  do  so.  Eventually 
all  the  A  ions  will  be  able  to  jump,  but  the  number  able  to  jump  at  any  one 
time  will  be  rather  small.  In  this  sense  the  number  of  defects  is  rather 
small— much. smaller  than  the  number  of  "vacancies"  in  the  crystal. 

Figure  5c,  also  after  Magneli,  shows  how  phases  of  general  formula 
A^Bg^  ^  can  be  built  up.  This  is  done  very  simply  by  allowing  the  zig-zag 

lines  of  extra  atoms  to  line  up  in  an  ordered  way  throughout  the  crystal. 

,  ./ 

The  phase  shown  is  A^B^  and  is  a  stoichiometric  compound  that  (as  drawn) 
contains  no  defects  at  all.  The  real  three  dimensional  phases  mentioned 
earlier  are  no  doubt  much  more  complicated  than  the  two  dimensional  case  con¬ 
sidered  here;  in  particular  both  cations  and  anions  will  be  involved  in  the 
ordering  process. 

To  summarize  the  discussion  so  far'  In  a  crystal  AB  having  a  structure 

n-x 

(as  revealed  by  x-ray  diffraction)  similar  to  AB  ,  the  stoichiometry  may  not 

n' 

be  a  reliable  guide  to  the  number  of  defects  able  to  diffuse.  For  example, 
a  vacancy  should  be  defined  as  a  site  such  that  if  on  atom  Jumps  into  t':at 
site  the  energy  of  the  crystal  is  unchanged,  m  this  sense  the  number  of 


vacancies  nay  be  much  less  than  the  stoichiometric  deficiency. 

The  example  we  considered  in  Figure  3a  is  still  a  non- stoichiometric 
crystal  (in  contrast  to  that  shown  in  Figure  3c)*  It  is  veil  known,  however, 
that  ve  may  have  a  stoichiometric  compound  with  the  same  sort  of  require¬ 
ments,  namely  that  a  given  ion  has  a  definite  number  of  nearest  neighbors, 
but  in  which  there  is  no  long  range  order.  A  very  well  known  example  is  that 

(26) 

of  ice.  Newman'  '  has  successfully  interpreted  the  structure  of  Oa0S^  and 

related  compounds  from  this  point  of  view.  Superficially  these  compounds 

have  a  defective  ZnS  structure,  but  Newman  showed  that  each  ion  in  the  crystal 

has  a  definite  coordination,  just  as  in  ice  each  oxygen  has  two  hydrogen  near- 

(27) 

est  neighbors.  As  Kroger'  '  has  pointed  out,  only  deviations  from  the  'ideal1 
short-range  order  in  which  a  given  ion  has  a  different  number  of  nearest 
neighbors  should  be  considered  a  defect  in  the  Frenkel  and  Schottky  sense. 

This  might  well  be  the  case  in  TiO^2^  and  VO^2*^  for  which  at  the  stoi¬ 
chiometric  composition  the  formula  based  on  a  rock  salt  structure  is  M0‘85°0’85* 
A  possible  explanation  is  that  each  metal  and  each  oxygen  has  5  nearest  neigh¬ 
bors  (rather  than  six  as  in  the  perfect  rock  salt  structure).  Migration  of 
(6ay)  oxygen  to  a  neighboring  site  requires  that  one  metal  will  have  four 
oxygen  neighbors  and  one  metal  six  neighbors.  These  'defects'  can  then  mi¬ 
grate  through  the  crystal. 

The  only  relevant  diffusion  work  known  to  the  author  is  that  of  Condit 

(50) 

and  "irchenall  on  Fe^S.'^'  These  authors  found  a  non-linear  dependence 
of  the  iron  diffusion  coefficient  on  the  degree  of  non-stoichiometry,  6;  vis. 

-  0.16  4  10"2  exp  (-U/rap  exp  {  -  (19  +  206)/Rt]>  (l8) 

« 

Although  ordering  of  defects  at  a  composition  near  Fe^Sg  reduced  the  iron 
diffusion  coefficient,  the  effect  was  less  than  one  order  of  magnitude 
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change  In  Dpe>  It  appears  to  the  present  author  that  this  is  strong  evi¬ 
dence  for  short-range  order  above  the  disordering  temperature.  The  compli¬ 
cated  dependence  of  L  on  6  (which  might  naively  be  equated  to  iron  vacancy 
concentration)  defeats  attempts  at  a  simple  explanation  based  on  a  random 
distribution  of  vacancies. 

Finally  under  the  heading  of  t;pe  0  materials  we  should  perhaps  include 

compounds  such  as  CUgS  and  AggS  in  which  an  almost  random  distribution  of 

ions  over  a  number  of  lattice  sites  occurs,  although  the  stoichiometry  may 

be  nearly  perfect.  In  CUgS  the  copper  atoms  are  in  three  sites  occupied 

about  8756,  35#  and  17#  respectively.^^  This  is  to  be  compared  with  a 

completely  random  distribution  in  which  33#  of  each  site  would  be  occupied. 

to  the  high  temperature  form  of  AggS  the  silver  is  thought  to  he  randomly 
(%2)  (tt) 

distributed.  '  Alien  and  Moore'  found  an  extraordinarily  high  diffusion 
coefficient  for  silver  in  this  phase,  viz. 

-  2.8  x  lO*1*  exp  (-3^50  cal/KT)  (19) 

AS*  calculated  from  equation  (3),  and  using  a  correlation  factor  J  =  0.50,^^ 
is  very  close  to  zero.  This  is  consistent  with  a  model  in  which  silver  is 
able  to  move  without  creating  any  additional  disorder.  Allen  and  Moore  com¬ 
ment  further  that  the  is  comparable  with  that  calculated  for  a  gas  of 
silver  atoms  at  the  same  temperature  and  density.  Hers,  then,  is  one  case 
where  complete  randomness  is  likely. 

Oxygen  Diffusion 

to  nearly  all  oxides  except  those  of  formula  MOg  (and  possibly  higher 
oxides)  oxygen  diffusion  is  considerably  less  rapid  than  cation  diffusion 
(R>0  la  one  exception).  This  situation  is  paralleled  by  that  obtaining  for 
the  halides  in  which  cation  mobility  is  normally  greater  in  the  lover  halides 
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but  anion  mobility  greater  in  (5^)  CdO  and  TiOg  it  appears  tbet 

anion  diffusion  occurs  by  a  vacancy  mechanism,  but  for  ZnO,  CUgO,  MgO,  NiO 

and  FbS  it  has  been  suggested  on  the  basis  of  the  pressure  dependence  of 

the  diffusion  coefficient  that  an  interstial  mechanism  is  operating.  An 

(yr) 

alternative  explanation  was  offered  by  Moorev-M/  and  more  recently  reviewed 
by  Woodbury.  This  theory  invoked  cooperation  of  cation  vacancies  in 
anion  diffusion  and  would  provide  an  explanation  of  the  observed  dependence 
of  D  on  a  positive  power  of  oxygen  pressure.  It  does  not  appear  possible 
to  distinguish  between  these  mechanisms  at  this  stage,  although  in  some 
Instances  a  general  indication  of  the  mechanism  of  oxygen  diffusion  can 
be  obtained  from  the  magnitude  of  the  pre- exponential  fuctor  in  the  dif¬ 
fusion  coefficient. 

Thus,  writing  D  *  Dq  exp  (  -  E^RT)  we  see  from  equation  (1),  (2)  and  (3) 
Do  -  (kT/h)  a 2  exp  {(AS*  *  AS*)/r} 

(kff/h)  e2  is  typically  about  10”  2  cm2  sec’1  at  1000“C,  and  ASi  is  gener¬ 
ally  expected  to  be  rather  small  for  diffusion  by  a  point  defect  mechanism/'*1' 
so  that  if  D0  i«*  very  different  from  10’2  cm2  sec”*  it  will  be  due  to  the 
factor  exp  ''vS'/R). 

For  experiments  at  constant  oxygen  (or  sulfur J  pressure  we  can  estimate 
the  order  of  magnitude  of  AS*  as  follows.  If  oxygen  vacancies  are  formed 
by  a  reaction  such  as 

2°0  — *  +  °2  (2°) 

AS  will  be  large  and  positive  including  not  only  the  contribution  from  the 
gaseous  oxygen  (6o,fc  cal  deg-1  mole’1  at  lJOO’K),  but  also  from  the  general 
loosening  of  the  lattice  in  the  proximity  of  the  vacancy.  This  might  con¬ 
tribute  a  factor  of  as  much  as  exp  (*»0  e.u/R  —  10^  to  the  pre- exponential 
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factor  D  .  On  the  other  hand  formation  of  interstitials  by  e.g. 
o 

2V-  +  0  - >  20t  (21) 

1  2  1 

will  be  expected  to  have  a  negative  AS  with  exp  (aS/r)  perhaps  ~  10 

Compare  for  example  (from  Table  2),  Dq  for  oxygen  diffusion  in  MgO  which 

(39) 

is  thought  to  occur  by  an  interctiai  mechanism'"^'  with  Dq  for  oxygen 

(to) 

diffusion  in  CdO,  which  almost  certainly  proceeds  by  a  vacancy  mechanism. v 

Finally,  in  connection  with  the  subject  of  oxygen  diffusion,  it  is  very 

interesting  to  note  that  recent  measurements  of  xenon  diffusion  in  several 

(4l) 

ceramic  oxides'  '  gave  results  surprisingly  similar  to  those  obtained  for 
oxygen.  While  this  may  be  fortuitous,  one  is  lead  to  speculate  on  possible 
mechanisms  of  diffusion  in  which  molecular  oxygen  (which  is  similar  in  size 
to  zenon)  plays  a  role.  Clearly,  more  experiments  of  this  type  would  be 
very  valuable. 

Interdiffusion  and  Reaction 

We  have  so  far  considered  only  self  diffusion  in  a  homogeneous  crystal 
as  measured  for  example  by  tracer  experiments.  However,  for  many  appli¬ 
cations  of  interest  to  the  chemist  and  ceramist  it  is  interdiffusion,  with 
possible  reaction  to  form  a  new  phase,  that  is  of  interest. 

Reactions  between  solid  oxides  to  form  new  phases  and  gas  solid  re¬ 
actions  are  the  subject  of  a  later  paper  so  that  they  will  be  omitted  here, 
the  present  discussion  being  limited  to  two  types  of  interdiffusion  process. 
Ihesc  are  (a)  the  interdiffusion* of  one  component  of  a  binary  compound  and 
the  compojnd  itself  (e.g.  Og  and  Cu„0)  resulting  in  a  change  of  stoichiometry 
end  (b)  the  interdiffusion  of  two  oxides  to  form  a  solid  solution  of  the 
same  structure  (e.g.  CoO  and  MgO). 

The  new  features  that  we  have  to  consider  are  that  now  ve  are  concerned 


-18- 

vith  diffusion  in  a  concentration  gradient  and  the  interdiffusion  process 
vill  involve  a  flux  of  at  least  two  different  species.  Until  the  advent  of 
isotope  tracer  techniques,  this  was  of  necessity  the  only  type  of  diffusion 
studied. 

Formally,  the  rate  of  interdiffusion  of  A  and  B  is  characterized  by  an 
interdif fusion  coefficient,  given  by'  ' 


3 


cadb  +  cbda 


(22) 


where  c^,  are  the  mole  fractions  of  A  and  B,  and  and  are  the 
diffusion  coefficients  of  A  and  B  in  the  solution  of  A  and  B.  It  may  be 
shown  tbat^2^ 


Dfi  *  D*  (l  +  dlnVg/ainCg) 

DA  c  Da  -1  +  dlnVd:LncA  ^ 


(23) 


Where  in  turn  D^,  Dg  are  the  diffusion  coefficients  of  A  and  B  as  measured 
in  a  homogeneous  solution  of  A  and  B  by,  for  example,  tracer  methods  (i.e. 
in  the  absence  of  a  chemical  concentration  gradient).  is  the  activity 
coefficient  of  A  in  the  solution. 

For  a  binary  solution  from  the  Gibbs -Duhem  equation  we  have  that 


ainVA  a  ainYp  a  SlnY 
^lnc'A  dlnCg  Bine 


(25) 


ao  that  from  (22)  we  have  the  Darken  Equation: 

IS 


cADg  +  cfiDA  (l  +  dlnY/dlnc ) 


(26) 


Strictly  equation  (26)  is  correct  to  the  extent  that  there  are  no 
correlation  effects  and  that  there  is  no  coupling  between  the  flow  of  A 
and  the  concentration  gradient  of  B  and  vice  versa.  For  a  more  complete 
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expression  for  D  reference  should  be  made  to  the  review  by  Howard  and  I/id- 


Turning  now  xc  some  specific  examples  let  us  suppose  that  B  is  a  gas 
that  can  dissolve  in  A  forming  a  dilute  solution.  If  A  is  considered  a  rigid 
matrix  then  and  as  the  solution  is  dilute  c^  ~  1  so  that  equation 

(23)  becomes 

B  «  D*  (1  +  ainY^ainCg) 


If  B  dissolves  in  A  as  ions  by  a  reaction  of  the  sort 


v-t 

B  — >  B  +  ve 


(27) 


Then  using  the  well-known  formalism  of  solution  theory 


(*3) 


V+l  V  V 

aB  -  a  «  C*Y+ 


(28) 


Vf 

where  Is.  the  mean  ionic  activity  and  the  +,-  subscripts  refer  to  B 

•  V+  • 

and  e  respectively.  Assuming  an  ideal  solution  of  B  and  e  we  have  Y+  =  Y„ 

and  vc,  «  c„  so  that 
+ 


(29) 


and  thus 


dlnag/djncg  «*  v  +  1 


‘5  «  (v  +  1)  D 


B 


(30) 

(3D 


As  an  illustration  we  may  take  the  diffusion  of  interstitial  zinc  in 

.  .  % 

zinc  oxide.  We  have  seen  that  zinc  dissolves  in  zinc  oxide  by  the  reaction 

“'g.s  — *  2»J  *  «’  UT) 

and  the  rate  of  uptake  of  zinc  can  be  obtained  from  conductivity  measure- 
(19) 

menta.  "  In  this  example  v  ■  1  so  that 

15  -  2  V*  (ZnJ) 


(32) 
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Note  that  D  (Zn£)  is  the  rate  of  diffusion  of  interstitial  zinc  in 
homogeneous  zinc  oxide  and  is  the  defect  diffusion  coefficient,  i.e,  D  (Znp  at 
D^/c  (Zn®)  (cf  equation  (b)). 

The  type  c '  diffusion  just  considered  is  more  frequently  discussed  from 

(kb) 

a  slightly  different  point  of  view.'  ‘  As  the  electron  is  more  mobile  than 
the  zinc  ion  it  will  tend  to  run  ahead  creating  a  ’diffusion  potential’  $ 
that  enhances  the  rate  of  diffusion  (accounting  for  the  factor  of  2  in  equa¬ 
tion  (32)).  The  flux  of  zinc  ions  of  charge  q  is  then  given  by 


J  (Zn®)  =  -  D*  {Zn®)  &c(Zn*)/dx  +  q  d§/dx  (33) 

and  the  flux  of  electrons  is  similarly 

J  (e” )  =  -  D*  (e")  dc(e”)/dx  -  q  d§/dx  (3*0 

Electrical  neutrality  requires  that 

J  (e“)  «=  J  (ZnJ)  (35) 

which  together  with  the  inequality 

D*  (Znz)«D*  (e~)  (3 6) 

and  D  =  D  (Zn®)  =  -J  (Zn®)/(dc(Zn®)/dx) 


yields  equation  (32).  This  type  of  diffusion  is emetimes  called  ambipolar 
diffusion. 

As  a  second  example,  we  may  consider  the  uptake  of  oxygen  by  a  crystal 

MO.  We  suppose  that  oxygen  is  taken  up  by  the  reaction  discussed  earlier 

.  •  * 

l/2  02  — >  0o  +  Vy  +  h®  (5) 

The  mobile  ionic  species  will  now  be  the  cation  vacancy  which  we  have 
assumed  ionized  os  shown.  In  thir.  instance  we  have 

3  -  2  l£ 


(37) 
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where  now  Dy  is  the  vacancy  diffusion  coefficient  Vv 

An  experiment  of  this  type  could  readily  "be  performed  as  follows.  A 
crystal  of  MO  at  a  low  oxygen  pressure  could  he  suddenly  subjected  to  greatly 
enhanced  oxygen  pressure  {all  the  time  staying  within  the  field  of  existence 
of  MO,  of  course).  The  rate  of  uptake  of  oxygen  (measured  by  the  change  in 
gas  pressure  or  perhaps  gravimetrically)  will  then  yield  the  vacancy  dif- 
flision  coefficient  directly.  This  is  a  method  that  deserves  more  attention 
than  it  has  been  accorded  so  far. 

There  has  of  necessity  been  a  simplified  account  of  interdiffusion  for 
two  simple  examples.  For  a  discussion  of  interdiffusion  taking  into  account 
cation  and  anion  mobility  and  Schottky  equilibrium  see  Brebrick. ' 

We  turn  now  to  a  second  type  of  interdif fusion  -  that  of  two  solids 
forming  a  solid  solution.  For  simplicity  the  discussion  will  be  limited 
to  on  greatly  oversimplified  example  to  illustrate  the  type  of  behavior  to 
be  expected.  The  example  considered  is  that  of  two  simple  oxides  AO  and  BO 
that  have  been  placed  together  and  allowed  to  interdiffuse. 

After  a  diffusion  anneal  the  interdiffusion  coefficient  nay  be  obtained 
from  a  Matano^^  analysis  from 

o 

?  =  *d<:  (38) 


where  the  origin  of  x  (the  Matano  interface)  is  defined  so  that 


x  dc  *  o 


(59) 


?>  measured  in  this  way  does  not  by  Itself  give  any  information  about  the 

*  * 

magnitude  of  or  Dg. 

In  order  to  find  out  what  type  of  concentration-distance  curve  to  ex¬ 
pect  we  make  the  following  simplifying  assumptions. 
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It  is  assumed  that,  oxygen  diffusion  is  negligible  and  that  cation  dif¬ 
fusion  is  by  a  vacancy  mechanism  and  further  that  the  concentration  of  vacan¬ 
cies  in  pure  AO  is  negligible  compared  with  that  in  BO. 

We  have  now  that 


®A  ‘  “vA  “V 


"bvb  cV 


(to) 


where  By^,  B^  measure  the  rate  of  vacancy  diffusion  by  Jump  into  an  A  01  a 
B  site  respectively.  It  will  not  affect  our  general  conclusions  ii  ve  let 
By »  Byg  “By.  We  then  have  very  simply 

*  -  IV  Cy  <*U 


Cy  will  depend  strongly  on  the  concentration  of  BO  (Cfi)  in  the  solution; 
suppose 


•  n 

Cy  “  Cy  cB 


(42) 


where  n  is  a  rational  number  of  the  order  of  unity.  Thus  from  (4l)  and 
(42)  if  B®  is  the  diffusion  coefficient  of  B  in  pure  BO: 

?  -  “v 

The  ejected'  concentration  profile  for  *5  of  this  type  can  be  obtained 
by  numerical  solution  of  Fick's  equation. ^  ^  Figure  4  shows  the  results 
for  n  »  2,  1  and  l/2.  In  this  diagram  c£  is  plotted  against  the  dimension¬ 
less  parameter 

X  -  */(AD-t)l/2  (W) 


The  curves  of  Figure  4  ere  rather  interesting.  They  show,  for  example, 
that  the  concentration  of  B  will  fall  to  zero  at-  a  finite  distance  into  the 
AO  so  that  there  will  be  a  definite  rate  of  advance  of  the  boundary  of  the 
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6 olid  solution  into  AO.  This  rate  is 

dx/dt  =  X0  (Dj/t)1/2  (45) 

where  X  is  the  value  of  X  at  which  falls  to  zero, 
o  B 

Although  it  involves  seme  rather  drastic  approximations  the  above  anal¬ 
ysis  probably  described  fairly  well  the  general  behavior  of  systems  such  as 
CoO-MgO  and  FeO-MgO.  A  recent  study  of  the  FeO-MgO  system ^  did  indeed 

show  a  concentration  profile  remarkably  like  that  of  figure  (4)  for  n  «  2. 

The  OoO-MgO  system^ ^  likewise  shows  a  sharp  boundary  between  (Co,  Mg)  0 

and  pure  MgO,  but  it  should  be  emphasized  that  despite  statements  to  the 
(49) 

contrary'  y/  the  rate  of  advance  of  the  boundary  is  not  a  measure  of  the 
rate  of  diffusion  of  cobalt  in  MgO  but  (cf  equation  (45))  perhaps  more  nearly 
measures  the  rate  of  diffusion  of  cations  in  pure  CoO. 

A  complete  analysis  of  a  system  such  as  this  presents  a  rather  diffi¬ 
cult  problem,  that  has  not  yet  been  solved  to  the  author's  knowledge.  One 
difficulty  that  might  be  mentioned  is  that  except  for  the  case  of  n  «  1 
(vacancy  concentrations  linearly  dependent  on  concentration  of  B)  the  total 
number  of  cation  sites  is  not  constant.  As  a  consequence  it  will  in  gen¬ 
eral  be  necessary  for  the  solid  system  as  a  whole  either  to  gain  or  to  lose 
oxygen.  This  will  be  a  particularly  important  effect  for  systems  involving 
oxides  such  as  FeO  where  the  vacancy  concentration  is  very  large. 

Diffusion  Bata 

Finally  we  present  in  condensed  form  a  summary  of  the  data  available 

.  .  •> 

on  diffusion  in  oxides  and  sulfides.  Only  data  directly  measured  by  tracer 
or  similar  techniques  have  been  Included.  No  effort  has  been  made  to  assess 
the  reliability  of  any  measurements,  as  this  would  be  extremely  difficult 
to  do  objectively.  It  must  be  pointed  out,  however,  that  in  a  few  instances 
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vhere  measurements  have  been  made  independently  on  the  same  system,  widely 
differing  results  have  been  obtained.  This  is  no  doubt  due  partly  to  the 
problem  of  impurities  mentioned  earlier;  however,  some  particularly  trouble¬ 
some  anomalies  remain.  One  notorious  example  is  that  of  the  diffusion  of 
zinc  in  zinc  oxide.  Only  two  representative  values  are  included  in  the 
table.  For  a  survey  of  this  problem,  reference  may  be  made  to  the  paper  of 
Moore  and  Williams  and  the  ensuing  discussion^^. 

In  a  great  number  of  cases  diffusion  coefficients  are  measured  in  a 
constant  atmosphere  (frequently  air)  with  temperature  as  the  only  parameter. 
D  is  then  reported  as 

D  *»  Dq  exp  (  -  E/RT) 

It  is  emphasized  again  here  that  it  would  frequently  be  very  profitable  also 
to  make  measurements  in  which  the  partial  pressures  of  the  ambient  gases 
are  systematically  varied.  Where  this  has  been  done  it  is  noted  in  the 
table. 

One  might  inquire  as  to  whether  there  were  any  significant  trends  in 

the  diffusion  coefficients  reported.  Lindner^ has  made  the  interesting 

observation  that  for  cany  oxides  the  cation  diffusion  coefficient  at  the 

-7  2  1 

melting  point  is  within  an  order  of  magnitude  of  10  cm  sec  ,  but  this 

rule  is  by  no  means  universally  true. 

The  enormous  range  of  reported  values  for  P  and  E  calls  for  comment. 

o 

There  appears  to  be  no  direct  correlation  between  Pq  and  F,  although  as 
might  be  expected,  in  general  a  Jarge  Dq  is  accompanied  by  large  E.  There 
is  a  more  direct  correlation  between  Pq  and  E/Ta;  in  figure  5,  E/Tn  is 
plotted  against  log  Dq.  It  may  be  seen  that  there  appears  to  be  an  approx¬ 
imately  linear  relation  between  these  two  quantities  for  the  oxygen  data, 
and  a  less  direct  correlation  for  calion  diffusion. 
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These  observations  present  an  interesting  theoretical  challenge. 
Ubbelohde^*^  has  remarked  on  the  very  large  pre- exponential  factor  in 
many  rate  processes  in  solids  and  was  led  to  suggest  that  this  might  be 
related  to  a  persistence  of  activation  energy  or  'overshoot1  phenomenon. 

By  this  means  an  ion  (particularly  an  interstitial)  could  Jump  several 
atomic  distances  at  once.  Alternatively  one  might  envisage  that,  rather 
than  diffusion  by  a  point  defect  mechanism,  there  was  a  cooperative  effect 
in  which  either  a  cluster  or  a  plane  in  the  crystal  moved.  For  a  recent 
discussion  of  the  relationship  between  activation  energy  and  pre- exponential 
factor  and  for  further  references  see  reference  (5*0. 
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TABLE  2 


Solid 

Diffusant 

_  2  -1 

D  cm  sec 

0 

^kcal 

Comment 

Reference 

cv 

Cu 

4.4  x  10-2 

36.1 

55 

0.12 

37.8 

56 

0 

6.5  x  Kf5 

39.3 

D  ~  P > 

°2 

57 

Ag 

6  x  10'5 

27.6 

58 

Zn 

4  x  10“4 

30.0 

59 

MgO 

Mg 

.21# 

79 

(60) 

0 

2.5  x  10‘6 

62.4 

\  ~  Po2  ^ 

39,  61 

Fe 

8.83  x  10“ 5 

1*1.7 

62 

Co 

5.78  x  10"5 

1*7.4 

62 

Ni 

1.80  x  10-5 

1*8.3 

62 

Zn 

1.1*8  x  10“5 

1*2.5 

63 

Xe 

~3x  10”6 

~  64 

4l 

NiO 

Ni 

2  x  10“4 

39.6 

64 

-2 

1.7  xio" 

56 

65 

0 

6.2  x  10"4 

57.5 

D  -  T  V2 

0  Og 

66 

ZnO 

Zn 

~  10 

73 

67 

1.3  x  10”5 

43.5 

50 

0 

6.5  x  1011 

165 

50 

In 

2,5  x  102 

73 

68 

FeO 

Fe 

1.1*  x  10’2 

30.2 

Fc01.08t 

69,  13 

CoO 

Co 

2.15  X  10’5 

34.5 

»o  ~  %  °° 

13 

FbO 

Fb 

105 

,  « 

66.6 

70 

0 

5-39 

22.4 

71 

CdO 

0 

8  x  106 

93 

D  ~  F 

°2 

4o 

BeO 

0 

68 

72 

Xe 

~  3  x  10”6 

~  64 

4l 

TABLE  2  (Continued) 


Solid 

Diffusant 

2  *1 

D  cm  sec 
o 

®kcal 

Comment  ReJ 

'erenco 

CaO 

Ca 

0.4 

81 

75 

BaO 

Ba 

-  IO30 

~  275 

above  1350*C 

74 

*2°3 

0 

1.9  X  IO3 

152 

T  >  l400#C 

75 

6.3  x  10"8 

57.6 

T  <  i4oo#c 

75 

A1 

28 

114 

76 

Xe 

3  x  10”6 

63.5 

4l 

Fe 

9.2  x  10”8 

27.6 

77 

er2°5 

Cr 

4  x  IO3 

100 

78 

0 

15-9 

101 

79 

.137 

61 

80 

Pe2°3 

Fe 

4  x  105 

112.0 

81 

1.3  x  106 

100 

82 

0 

10 P° 

146 

83 

Y2°3 

Y 

2.4  x  lo"4 

43.9 

84 

UCL 

U 

0.23  k 

104.6 

H,0  atmosphere 

85 

2 

4.3  x  10 

88.0 

H 2  atmosphere 

86 

others 

t— 

CO 

7X°2 

0 

5.1  x  10”3 

29.8 

(15*  CaO) 

10 

Xe 

6.5  x  10”8 

42 

4l 

Ti°2 

0 

2  x  iO"3 

60 

(.0156  AlgOj) 

88 

Fe 

2  x  10”2 

55-0 

in  air 

T7 

1.9  x  10’ 1 
.  « 

55.4 

in  vacuum 

77 

SnO 

2 

3n 

io6 

1 19 

78 

°€°2 

0 

io”7 

28.5 

89 

(10‘3) 

27.5 

90 

Fe 

5.2 

55.0 

F>2.993  °4 

69 

6  x  io5 

84 

r  -  0-^(9  0 

91 

TABLE  2  (Continued) 


Solid 

Diffusant 

2  -1 

D  cm  sec 

0 

h«l 

Comment 

Reference 

CaSiO* 

3 

Ca 

7.1*  x  104 

112.0 

81 

PbSiO, 

j 

Pb 

85 

59.5 

8l 

PbSiO^ 

Fb 

8.2 

47 

81 

SrTiO, 

3 

0 

~10'5 

~  27.5 

independent  of 

ZnPegO^ 

Zn 

8.8  x  10* 

86.0 

93 

Fe 

8.5  x  102 

82 

93 

ZnCr2°4 

Zn 

6c 

85.5 

78 

ZnCr204 

Cr 

8.5 

81.0 

78 

CoCr204 

Co 

10"5 

51 

94 

Cr 

2 

72.5 

94 

NiCr204 

Ki 

1.5  x  10'5 

6l.4 

52 

Cr 

74 

72.5 

78 

0 

.017 

65 

83 

ZnAlgO^ 

Zn 

2.5  x  10'2 

78.0 

78 

NiAlgO^ 

Ni 

2.9  x  10“5 

53.5 

52 

SnZn^O^ 

Sn 

2  x  105 

109 

78 

Zn 

37 

76.3 

78 

CaFe  C. 

2  4 

Ca 

30 

86 

81 

Fe 

0.4 

72 

61 

Wl2 

0 

0.4 

65.4 

95 

A^S 

Ag 

2.8  x  lo’4 

3.5 

in  cubic  AggS 

33 

S 

2.4  x  10'4 
,  » 

25 

n 

98 

FbS 

Fb 

8.6  x  10“5 

35 

•toichiometrie 

97 

FbS 

S 

6.8  x  lo"^ 

31.5 

II 

9L 

FeS 

Fe 

•ee  text  equation  (l8) 

30 

3 

10»-17 

240 

30 

TABLi.  2  (Continued) 


Solid 

Diffi’sanc 

_  2  -1 

D  cm  sec 

0 

\c&l 

Comment 

Reference 

NiS 

Ni 

1.1  x  10“2 

25-5 

|i  c-axis 

99 

8.5  x  io-5 

25.6 

X  c-axis 

s 

2.5  x  10S 

65.5 

fl  c-axis 

99 

2.2  x  106 

87.5 

X  c-axis 

CdS 

Cd 

3*^ 

k6 

in  6at.  Cd  Vapor 

38 
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Introduction 

A  large  number  of  chemical  reactions  of  considerable  industrial  importance 
involve  the  diffusion  controlled  transport  of  matter  in  oxides.  As  examples  we 
may  cite  the  oxidation  of  metals,  the  reduction  of  oxides,  sintering  of  ceramics 
and  the  formation  of  ternary  oxides  (such  as  spinels)  from  binary  oxides.  It  is 
clearly  of  considerable  interest  to  be  able  to  relate  the  diffusion  coefficients 
characterizing  these  processes  to  the  diffusion  coefficients  obtained  from  radio¬ 
isotope  studies  in  homogeneous  materials. 

In  order  to  clarify  the  definitions  used  and  also  to  make  clear  the  under¬ 
lying  assumptions,  we  first  sketch  a  derivation  of  the  relationship  between  the 
'chemical1  and  'tracer1  diffusion  coefficients,  following  the  general  treatment 
given  for  example  by  Darken.1 

In  en  experiment  in  which  a  vet  flux  of  matter  is  observed,  a  chemical  dif¬ 
fusion  coefficient  can  be  defined  in  terms  of  the  experimentally  measured  quan¬ 
tities  by  Pick' 8  equation 

\  =  -J^/grad  ci  (1) 

In  practice  an  integrated  form  of  Fick's  'second*  equation  is  usually  used  but 
this  is  equivalent  the  definition  (l), 

To  relate  D  to  other  measuiabl*-  quantities  it  is  assumed  that  the  driving 
force  for  isothermal  diffusion  is  the  gradient  of  electrochemical  potential  p, 
so  that  the  velocity  of  the  diffusing  species  I  is 

Vi  *  Jj/Gi  *  Srad  Hj.  (2* 

where  u^  is  the  mobil'ty  of  species  i.  Wilting 

o  _  , 

■  u^  ♦  RT  Ir. 


(3) 
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we  have  from  eqs.  (l)  and  (2) 

^  =  UA  HT  (din  a,/dln  ct)  (4) 

In  order  to  relate  D^  to  the  tracer  diffusion  coefficient,  D^,  It  is 
supposed  that  in  a  dilute  solution  the  tracer  atom  behaves  identically  to  the 
solvent  atoms  and  that  din  a^/dln  c^  is  unity.  However  it  is  possible  that 
atomic  jumps  are  not  completely  random  but  correlated  so  that  we  write 

D*  =  f  ^  RT  (5) 

where  f  it  correlation  factor.  We  return  to  a  discussion  of  f  below.  Sub¬ 
stitution  of  eq.  (5)  into  eq.  (4)  yields  finally  the  Darken  equation 

Uj  =  (l/f)  D*  (din  aj/dln  ci)  (6) 

Wlistite  (FeO  )  would  appear  to  be  an  ideal  material  with  which  to  test 

p 

eq.  (6).  It  has  been  known  since  the  experiments  of  Pfeil  that  oxygen  dif¬ 
fusion  in  negligibly  small  in  this  phase.  Several  concordant  measurements  of 

tracer  diffusion  have  been  made"  the  results  of  some  of  these  measurements 

5  5  * 

are  hown  in  Figure  1.  It  has  also  been  shown'"  that  D^e  increases  with  in¬ 
creasing  x  in  FeOx  especially  at  higher  temperatures  (above  9C0°C);  this  is  a 
significant  observation  to  which  we  return  below. 

6  7  8 

There  have  also  been  a  number  of  measurements' of  activity  as  a  function 
of  concentration  which  are  again  in  good  agreement  with  each  other.  From  these 
we  have  calculated  that  at  1050®C  in  FeOx 

din  aFe/dln  cFe  -  28.5  x2  (7) 

Levin  and  Wagner^  studied  chemical  diffusion  in  this  system  by  partial 
reduction  of  FeOx.  7keir  results  were  in  fair  accord  with  eq.  (6)  but  interest¬ 
ingly  it  was  found  that  *&  decreased  with  increasing  x  whereas,  If  f  is  constant, 
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eq.  (6)  predicts  an  increase  of  ?>  with  x  (both  and  din  a^dln  Cpe  increase 
with  increasing  x). 

We  decided  to  examine  this  system  a  little  more  closely,  to  see  whether  we 
could  resolve  this  discrepancy;  as  will  be  seen  our  results  are  in  even  less 
accord  with  eq.  (6). 

Experimental 

FeOx  was  made  in  situ  in  a  vacuum  microbalance  apparatus  by  oxidizing  1  mm 
thick  99.999$  (with  respect  to  metallic  impurities)  iron  from  United  Mineral 
and  Chemical  Corp.  in  an  approximately  equimolar  CO^/CO  mixture.  The  material 
was  always  maintained  within  the  field  of  stability  of  wttsuite  after  oxidation. 
Examination  of  the  specimen  after  a  series  of  experiments  showed  that  it  was 
coarsely  crystalline  with  grains  several  mm  in  size  and  with  a  hollow  core  result¬ 
ing  from  the  outward  diffusion  of  iron  during  oxidation.  After  equilibration 
in  a  given  CO^/CO  mixture  the  pressure  was  reduced,  to  less  than  1  mm  (oxygen 
activity  being  unchanged  during  this  process).  A  CO^/CO  mixture  of  slightly 
different  oxygen  activity  was  then  added  and  the  kinetics  of  oxygen  uptake  or 

loss  by  the  sample  was  followed.  The  diffusion  coefficient  was  calculated 

0 

utilizing  the  well-known  solution  for  this  diffusion  problem  from  the  data 
for  the  last  30$  of  the  weight  change.  The  total  weight  change  was  typically 
several  tenths  of  a  milligram  and  could  be  measured  with  a  precision  exceeding 
1$.  It  is  characteristic  of  this  type  of  measurement  that  only  relative  weight 

changes  need  be  measured  but  in  fact  the  abcolute  weight  changes  agreed  well 

6  7  8 

with  the  published  data.  9  9  Thp  thermocouple  was  calibrated  in  situ  by 
observing  the  melting  point  01  a  sample  of  pure  gold  hanging  in  place  of  the 
Supple  • 
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Results 

Most  of  the  measurements  were  confined  to  constant  temperature  of  1055°C. 
Results  for  two  series  of  measurements  made  on  two  different  samples  are  shown 
in  Figure  2.  The  general  behavior  of  I$pe  decreasing  with  x  as  reported  by  Levin 
and  Wagner  (and  in  subsequent  measurements  reported  by  Prof.  Wagner  at  this 
symposium)  are  confirmed  but  our  diffusion  coefficients  are  consistently  lower 
by  about  an  order  of  magnitude.  A  possible  explanation  of  this  difference  is 
given  below. 


Discussion 

Representative  result 3  presented  in  Table  1  illustrate  clearly  that  the 
Darken  equation  with  f  »  1  is  not  obeyed  in  this  system.  We  feel  that  there 
are  two  possible  explanations  for  this  result. 

1.  When  the  composition  of  FeO^  is  suddenly  changed  near  the  surface  there 
is  a  change  in  lattice  parameter.  The  shape  of  the  sample  will  not  be  able  to 
change  to  accommodate  this  change  in  molar  volume  so  that  there  will  be  a  stress 
gradient  across  the  sample.  To  relieve  this  stress  there  will  have  to  be  a 
flux  of  matter  normal  to  the  concentration  gradient.  If  the  stress  gradient  is 
sufficiently  large  it  is  entirely  possible  that  the  rate  of  change  of  composition 
is  controlled  by  the  rate  of  relaxation  of  stress,  i.e.  by  diffusion  of  iron  in 
the  e.bsence  of  a  concentration  gradient  (it  is  assumed  that  oxygen  can  be  trans- 

4 

ported  in  the  gas  phase)  as  for  example,  in  sintering. 

In  terms  of  the  equations  given  above  the  chemical  potential  appearing  in 


eq.  (5)  must  be  replaced  by 


10 


*Fe  “  ^Fe  +  131  10  ‘Fe  + 


(8) 


where  a^e  is  the  activity  of  iron  in  the  unst'^ssed  material  and  is  the 
change  in  chemical  potential  due  to  the  stress.  Explicit  expressions  for 
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ln  terms  of  the  stresses  and  elastic  coefficients  have  been  given  by  Darken,  Li 
and  Oriani.'5'0  Equation  (6)  becomes  now 


-a 

5.  -  B,  K  «1  l  a(6iJi) 

1  Y  din  c1  Iff  din  c± 

l 


Only  the  first  term  in  square  brackets  is  obtainable  from  equilibrium 
measurements.  Evaluation  of  the  second  term  requires  a  detailed  knowledge  of 
the  elastic  compliances  of  the  crystal  but  it  is  important  to  note  that  it  is 
of  the  opposite  sign  to  the  first  term.  We  may  illustrate  this  last  point  by 
a  veay  rough  approximation  (not  to  be  taken  too  seriously,  but  simply  given  to 
show  that  the  effect  of  stresses  must  be  taker,  into  account)  as  follows.  If  the 
shear  modulus  is  much  greater  than  the  bulk  modulus  we  may  suppose  the  stress 
to  be  essentially  a  hydrostatic  pressure  P  given  approximately  by 
P  -a  "(l/p)(dV/v)  where  f)  is  the  compressibility  and  AV/V  the  fractional  change 
in  volume  of  the  material  that  would  occur  on  changing  the  concentration.  6|j 
is  then  -(V1/0)(av/v)  and  finally 

~Yi  auv/v)  (10) 

RT  din  ci  8HT  din  c1 

11  12 

From  the  data  of  Jette  and  Foote  and  Levin  and  Wagner  relating  lattice 
parameter  to  composition  we  estimate  VFe  »  5.0  cm^  mole-1,  d(AV/v)/dln  cFe  =  O.k, 
so  that  with  0  ~  10~^  atm-1  the  right  hand  side  of  eq.  (10)  is  equal  to  about 
-20  compared  with  about  +30  for  din  &yj din 

This  consideration  is  particularly  important  for  a  material  of  wide  compo¬ 
sition  ra^ge  such  as  wiistite.  For  an  oxide  like  NiO^  for  example  din  a/dln  c 
is  of  the  order  of  10^  30  that  the  second  term  in  eq.  (9)  will  t :  relatively  un¬ 
important  . 


-6- 


2.  We  have  still  to  explain  the  observation  that  decreases  as  x  in¬ 
creases.  This  can  be  explained  if  f  Increases  with  increasing  x,  as  is  quite 
possible  in  a  \^ghly  defective  material  such  as  vttstite. 

In  a  system  containing  a  dilute  solution  of  vacancies  Jumping  strictly  at 
random  f  will  be  less  than  one  due  to  the  fact  that  tracer  Jumps  are  not  random 
but  sequentially  correlated  with  a  higher  tnan  statistical  probability  of  an 
immediate  return  Jump.  If  however  one  has  a  concentrated  system  of  defects 
their  Jumps  will  no  longer  be  random  but  on  the  contrary  it  is  likely  that  there 
is  correlation  between  two  (or  more)  simultaneous  atomic  Jumps.  This  will  be 
because  of  interactions  between  defects  such  as  vacancies  making  certain  con¬ 
figurations  more  stable.  Simultaneous  Jumps  that  preserve  the  integrity  of  the 

* 

complex  will  then  be  preferred  over  randan  Jumps. 

It  is  quite  likely  then  that  the  decreasing  1)  observed  in  this  system 
reflects  an  increase  of  f  from  O.'fQ  for  a  dilute  solution  of  vacancies  in  the 
(hypothetical)  material  with  x  close  to  unity  to  a  value  somewhat  greater  than 
one  in  material  with  x  ~  1.1.  It  is  possibly  significant  in  this  connection 
that  D  appears  to  level  off  at  a  composition  close  to  the  W1  -  W2  transition 

7 

of  Baccah  and  Valiet  which  is  likely  to  represent  an  ordering  of  defect  clusters. 
Suranary 

Chemical  diffusion  in  the  FeOx  (wtistite)  system  indicates  that  the  Darken 
equation  does  not  correctly  predict  the  relationship  between  chemical  and  tracer 

# 

It  is  possibly  helpful  to  consider  tv.  extreme  examples .  m  a  direct 
place  exchange  mechanism  of  diffusion  by  tracer  atoms  a  net  flux  of  matter 
(or  ionic  conductivity)  is  not  possible  so  that  f  -  ®.  vay  of  contrast  the 
’caterpillar'  mechanism  proposed  by  Yokota^  and  vhicn  consists  essentially  of 
simultaneous  linear  Jumps.  The  net  flux  of  atoms  can  be  much  greater  than  that 
of  tracers  so  that  f  «  1  (Note  that  f  defined  by  Yckota  is  the  reciprocal  of 
that  used  in  this  paper). 
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diffusian  coefficients.  A  modification  of  the  Barken  equation  is  proposed 
which  takes  into  account  the  stresses  present  in  an  inhomogeneous  material. 

The  results  also  indicate  that  in  very  defective  materials  (such  as  vtistite) 
the  correlation  coefficient  may  veil  be  concentration  dependent. 

Research  sponsored  by  the  Air  Force  Office  of  Scientific  Research,  Office 
of  Aerospace  Research,  under  grant  NrAF-AI’OSR-719-65* 


Table  1 

Observed  and  calculated*  chemical  diffusion  coefficients  in  FeO 

x 


T°C 

X 

lo6  D  (calc) 

’O^  D  (obs) 

Author 

1100 

1.125-1.05 

2.7 

3.2 

Levin  and  Wagner 

1100 

1.125  1.10 

6.7 

0.87 

II  IT  It 

1055 

1.07 

2.5 

0.15 

This  Work 

1055 

1.13 

6.1 

0.06 

»  »t 

* 

Calculated  from  eq.  (6)  with  f  •*  1. 


log  D 


FeO 


1.09 


-7r  *  \. 


- —  Dfe  (cm2 sec'1) 

=  (4.M0'3)exp(-27.8  kcal/RT) 


*  Carter  &  Richardson 

o  Desmarescaux  &  \Lacombe 

*  Himmel,  Mehl  &  Birchenall 

*  Pettit  (oxidation) 
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Diterdiffusian  Coefficients  from  Conductivity  Measurements 


-  Application  to  CUgO 

R.  H.  Campbell,  W.  J.  Kass  and  M.  O'Keeffe 
Chemistry  Department,  Arizona  State  University,  Tenrpe,  Arizona 


Introduction 


In  many  non- stoichiometric  oxides  with  Ionized  vacancies,  the  conductivity 

\ 

is  proportional  to  the  departure  from  the  idee  L  composition.  This  is  the  case 
for  example  in  CoO.1  In  CUgO  the  cond. ctivity  is  accurately  proportional  to 
the  square  root  of  the  vacancy  concentration  (as  shewn  by  the  oxygen  pressure 
dependence  of  conductivity  and  vaeancy  concentration').  It  would  appear 
possible  therefore  to  use  conductivity  as  a  measure  of  composition  of  the 

li 

sample,  and  indeed  several  authors  have  used  this  fact  to  follow  the  kinetics 
of  interdiffusion  in  oxides. 

It  is  necessary  however  to  proceed  with  caution  as  in  general  the  apparent 
conductance  of  an  inhomogeneous  sample  is  not  proportional  to  the  average  con* 
centration  of  defects.  Provided  that  there  is  local  electrical  neutrality, 

(i.e.  in  the  absence  of  a  space  charge)  it  is  possible  however,  to  relate  the 
local  conductivity  a  to  the  local  concentration  of  defects,  c. 

Theory 

We  have  considered  the  case  of  a  rectangular  prism  "ai  <  ^ 

“a2  <  x2  <  a2*  “a3  <  x3  <  a3  Poroue  electrodes  on  the  x^,  faces  so 
that  the  conductance  can  be  measured  in  the  x^  direction.  We  suppose  that  the 
initial  conductivity  is  due  to  an  initial  homogeneous  concentration  of 
defects  For  times  greater  than  zero  the  surface  concentration  is  held  at 
cf  so  that  the  conductivity  comes  tc  a  value  of  as  t  -  The  concentration 
at  any  point  in  the  material  is  then  given  by  the  well-known  expression 

c  (x^  Xg,  x?)  -  cf  -  (cf-ci)  S  (5X)  S  ($2)  3  (?5)  (1) 

n 

S  (5l)  ■  (2/tt)  I  ^(-l)n/{n+l/2)J  exp  ["-B^n+l/S)2’''2]  cos  (2) 


where 


2- 


and  In  turn 


2 

Si  *  x±/*±  >  Pi  ■  Dt/e^ 


The  local  conductivity  is  a(x^f  x2>  x^)  and  as  discussed  above  ve  have 
a(x]_,  Xg,  x?)  =  a[c(xx,  Xg,  XjjJ®  (k) 


with  A  a  constant  and  m  usually  1  or  l/2. 

The  apparent  conductivity,  o,  measured  in  the  x^  direction  is  obtained  by 
multiplying  the  conductance  of  the  sample  by  the  length/area  ratio  and  is  re¬ 


lated  to  the  local  conductivity  by 


a  "  2  a. 


-  f- 

ia2  4,  r& i 


5  L  L  a  ax*  a* 


"bl  a2 


2  “*1 


From  eqs.  (5),  (4)  and  (l)  one  obtains  after  seme  manipulation 


where  I  =  £  [l  +  -l)  S  (^J  S  (|g)  S  (§5)]“  §x  d§2 

end  a  •  o^/c^, 

It  is  of  interest  to  compare  2}  as  given  by  eq,  (6)  with  the  corresponding 
quantity  W  ■  (wj-w)/(wf.-vi)  where  vf  and  are  the  final  and  initial  amounts 
of  defects.  The  explicit  expression  for  W  Is  given  in  reference  h. 

In  Figure  1  ve  show  W  and  I  as  a  function  of  the  dimensionless  quantity 
Dt/a|  for  various  shape  samples  and  values  of  o^/c^  for  the  case  m  *  1.  It 
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may  be  seen  that  in  general  W  and  £  are  very  close  to  each  other  for  long,  thin 
samples  when  o^/a^  is  not  very  different  from  unity.  In  other  cases  hew  ever 
the  assumption  that  £  ■  W  may  lead  to  large  errors  in  the  calculated  diffusion 
coefficients . 

Application  to  Cuprous  Oxide 

Attempts  to  study  interdiffusion  in  cuprous  oxide  using  a  vacuum  micro- 
balance  technique  were  severely  limited  by  the  nature  of  the  material.  Spe¬ 
cifically,  at  temperatures  above  1000“C  evaporation  is  very  rapid  so  that  the 
small  weight  changes  upon  change  of  stoichiometry  are  hard  to  determine.  At 
lower  temperatures  the  cuprite  phase  field  is  very  narrow,  again  precluding 
reliable  measurements.  For  these  reasons  we  decided  to  resort  to  electrical 
measurements  in  this  system. 

Two  single  crystal  specimens  were  used;  sample  1  was  4.50  x  4.65  x  6.34  mm 
and  sample  2  was  I.91  x  1.91  x  15*9  mm.  lb  each  case  the  conductivity  was 
measured  along  the  longest  direction  using  two  probes.  After  equilibration  at 
low  oxygen  pressure,  oxygen  was  rapidly  admitted  to  the  sample  and  the  con¬ 
ductivity  monitored  as  a  function  of  time.  Strictly,  in  order  to  use  the 
expressions  derived  above  for  E,  porous  contacts  covering  the  whole  face  of  the 
crystal  should  be  used.  The  use  of  smaller  contacts  is  likely  to  be  less  serious 
in  the  case  of  sample  2.  As  conductivity  changes  are  very  rapid,  the  sample 
current  with  known  constant  applied  voltage  was  recorded  with  a  recorder  speed 
of  6  inches/minute.  Plots  of  E  versus  log  t  were  superimposed  on  theoretical 
plots  of  E  versus  log  (Dt/a^  )  calculated  according  to  eq.  (6).  '[he  difference 

A 

in  the  abscissa  of  the  two  plots  is  equal  to  log  (E/a*)  so  that  D  is  readily 
calculated.  As  the  conductivity  is  very  temperature  sensitive  it  is  quite 
possible  to  measure  erroneous  values  of  however,  this  becomes  evident  when 
the  experimental  curve  is  coopered  to  the  theoretical  curve,  and  any  experimental 
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rasults  differing  significantly  from  the  theory  were  rejected. 

Because  of  the  very  rapid  diffusion  observed  there  was  considerable  scatter 
of  the  experimental  results.  These  are  reported  in  Table  1  with  results  for 
each  crystal  at  each  temperature  reported  as  an  overage.  The  results  for  sample 
2  were  in  general  better  and  are  in  fair  agreement  with  the  theoretical  values 
(discussed  below).  The  results  for  sample  1  are  too  high;  it  is  not  clear  at 
present  whether  this  is  due  to  the  fact  that  the  experimental  conditions  did 
not  conform  to  the  boundary  conditions  used  in  deriving  eq.  (6)  or  whether 
effects  such  as  surface  diffusion  are  significant. 

Theory  for  Cuprous  Oxide  and  Comparison  with  Experiment 

We  can  calculate  the  expected  interdiffusion  coefficients  as  follows. 

5 

Neglecting  oxygen  diffusion,  the  rate  of  interdiffusion  of  oxygen  and  CUgQj^ 
will  be  given  by  the  Darken  equation 

D  =  din  a/ils  c  (7) 

where  is  the  tracer  diffusion  coefficient  of  copper  in  cuprous  oxide. 
According  to  Moore  and  Selikson^  at  an  oxygen  pressure  of  0.1  torr 

iogio  -  -  1.36  -  79CO/T  (8) 

Under  the  same  conditions  it  was  found  by  O'Keeffe  and  Moore^  by  direct 
gravimetry  that 

log1Q  x  -  0.02  -  U760/T  (9) 


* 

Objections  raised  in  the  previous  paper  to  the  Darken  equation  are  not 
valid  here.  This  is  because  now  din  a/dln  c  is  very  large  (~  10^)  and  we  are 
considering  a  dilute  solution  of  vacancies  (<  0.l£). 


5- 


9  In  x/3  In  P^  -  0.27  (10) 

Nov  setting  din  a/dln  c  ■  l/2  din  P^dln  (l+x)  and  recalling  that  x  «  1 
ve  obtain  from  eqs.  (7)  -  (10) 


log^  d  -  -1.U  -  jito/T  (n) 

'fable  1  also  records  15  as  calculated  from  eq.  (11).  The  agreement  with 
sample  2  Is  really  quite  remarkable  when  it  Is  remembered  that  the  three  separate 
measurements  being  compared  are  all  subject  to  large  errors.  Perhaps  the  most 
difficult  measurement  is  that  of  the  non- stoichiometry  x  and  it  is  gratifying 
that  the  interdif fusion  results  do  confirm  the  order  of  magnitude  at  least  of 
the  previously  reported  results. 

The  15  reported  here  may  be  interpreted  as  a  vacancy  diffusion  coefficient. 

The  very  high  diffusion  coefficients  are  in  accord  with  the  observation  that 

internal  equilibrium  can  readily  be  reached  in  cuprous  oxide  down  to  very  low 
2  7 

*  It  is  worth  remarking  also  that  in  general  cation  mobility 


temperatures . 


is  very  nigh  in  the  chalcogenides  of  the  univalent  metals. 

This  research  was  sponsored  by  the  Air  Force  Office  of  Scientific  Research, 
Office  of  Aerospace  Research,  United  States  Air  Force,  under  AFOSR  grant 
Nr.  AF-AF0SR-717-65. 


Table  1 

Measured  and  Calculated  Values  of  log^Q  75  in  CugO 


T°C 

* 

Sample  1 

Sample  2 

Calculated 

933  +  1 

-3.4  (1) 

-3.6  (2) 

-3-71 

973  i  1 

-3.15  (5) 

-3.4  (2) 

-3.63 

1001  +  1 

KN 

CVJ 

• 

K~\ 

• 

-3.75  (4) 

-3.57 

1032  +  1 

-2.6  (2) 

-3.45  (3) 

-3.51 

1064  +  l 

-2.5  (2) 

-3.4  (5) 

-3.46 

Numbers  in  parentheses  are  the  number  of  separate 
measurements  at  that  temperature. 


Calculated  from  eq.  (ll) 


Figure  1  -  Theoretical  curves  for  conductivity  and  weight  change  as  a 
function  of  time  for  various  shape  samples.  The  numbers  on 
the  curves  are  the  values  of  cr^/o^.  The  curves  labeled  W 
refers  to  weight  change,  the  others  to  conductivity  change. 
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A  Fortran  Program  to  calculate  conductivity  and  weight  as  a  function  of  tine 
during  lnterdlf fusion.* 

The  following  program  la  designed  to  calculate  £  and  V  as  described  In  the 
previous  paper.  The  necessary  Input  data  are  provided  on  three  cards, 
card  (a)  title  in  alphanumeric  form  In  columns  l-8o 
card  (b)  m  as  a  floating  point  number  In  columns  1-10 
card  (c)  the  values  of  a^/a^,  a^/a^  and  cr^/a^  in  column  1-10,  11-20, 
21-30.  For  the  most  useful  range  of  output  should  be  *  a^. 

This  set  of  three  cards  may  be  repeated  for  each  sample;  the  final  card  after 
all  the  sets  should  be  blank. 

The  output  Is  BETA?  a  a  Dt/a^  ,  log^  (B^),  £  and  W.  In  the  heading 
LZ/LX  ■  a^/a ji  LZ/LY  ®  «j/a2’ 

The  program  has  been  found  to  ran  well  for  m  «  0-5  and  1  (In  the  latter 
case  a  mac*  simpler  program  could  be  devised)  and  for  <  1.  The  program 

Is  doubtless  Inefficient  but  the  demand  In  computer  time  Is  modest  so  that 
little  effort  has  been  made  to  streamline  It.  A  series  of  calculations  for 
one  sample  takas  about  five  minutes  using  a  CDC  3^00  computer. 


*  Written  by  V.  J.  Kmo. 
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PROGRAM  PIFCON 

C  PROGRAM  TO  DETERMINE  diffusion  coefficients  from  CONDUCTIVITY 
C  MEASUREMENTS 

C  BRUTE  FCRCF  INTEGRATION 

C  OMEGA  MIST  BE  LESS  THAN  1  (GAS  UPTAKE) 

C  OMEGA  s  SIGMA  INITIAL  /  SIGMA, FINAL 

C  EM  =  DEPENDENCE  OF  SIGMA  ON  CONCENTRATION,  THAT  IS, 

C  SIGMA  s  CONSTAMT*CONCENTRATION **EM 

C  A3?  =  A3/A?  LENGTH  OF  Z  AXIS  OF  CRYSTAL  /  LENGTH  OF  LONGEST  0 

C  X  OR  Y  AXES 

C  A31  s  A3/A1  LENGTH  OF  Z  AXIS  OF  CRYSTAL  /  LENGTH  OF  SHORTEST  0 
C  X  OR  Y  AXES 

C  A31  SHOULD  BE  A  LARGER  NUMBER  THAN  A32 

COMMON  SK ( 50 ) »  PROD(50),  A31,  A32,  BET  A3 ,  XI3,  OMEGA.  EM 

DIMENSION  T I  TLE ( 20 ) 

LOGiOF(X)  =  L0GF(X)/2,;5Q25851 
PI  =  3.1415926536 

944  READ  970,  (TITLE(I),  I  =  1,  20) 

970  FORMAT(20A4) 

Read  900,  EM 
42  FORMAT(lHl) 

IF ( EM )  543,  942,  943 
943  READ  900,  A31,  A32,  OMEGA 
9C0  FORMAT(3Fin.O) 

PRINT  972 
972  FORMAT ( 1H1 ) 

PRINT  970. (TITLEC I),  I  =  1,  20) 

PRINT  oiO,  EM,  OMEGA 

9j0  FORMAT (  5H  M  =,  F6.3,10X,  19HSIGMA  I  /  SIGMA  F  =,  F7,4) 

PRINT  911.  A  3 1 ,  A32 

911  FORMAT  ( 9H  1.7/LX  =,  F8.4,  4X,  7HLZ/LY  =,  F8.4  //) 

PRINT  926 

926  F0RMAT(6X,  5  H  R  E  T  a  3 ,  5X,  9HL0G  (  BET  A  )  ,  5X,  |jHS  I  GM  A ,  6X,  A  H  R  E  I  G  H  T  ) 

A  =  10 . 0*LOGl HF( A31 ) 

KLOG  =  A 
AA  =  Kt.CG 
AALCG  =  A  A/ 10.0 
PETLOG  =  -4.5  -  2. 0*AA|  OG 
90?  DO  9.J0  MN  2  50 

PETLOG  2  HETLOG  ♦  0.1 
B  E  T  A  3  =  EXPHUETI. OG*2, 3025851) 

PROD ( 1 )  2  o.O 
CALL  SjMPSO\i<SjHP) 

SIOMA  2  U.n/ (1.0  -  OMT^A  )  )  *  (  1  .  0  -  (1.0/SlMP)) 

PE  T  A  1  ;  A  3  1.  ♦  A  3 1  * H E  T  A 3 
BETA?  ;  A3?* A ??*META3 
CALL  S  IVPW  r  (  META* .  i.'TX,  EH) 

I  F ( A  3 1  -  A3?)  960 ,  9  >1,  960 
96i  wry  =  w T x 
GO  To  7 

96Q  CAI.L  S  I  MPW  T  ( t*E  T  A  ?  ,  WTY,  R|) 

967  IK  A  <1  -  l.f!)  96?,  963,  96? 

063  HT/  -  yTX 
GO  TO  96S 

96?  IF ( A3?  -  1.6)  064,  9*6,  964 
9 6.6  WT/  1  W  T  Y 


< 1. rt/SlMP)  ) 


c 


r 
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GO  TO  965 

/  "  964  CALL  SIMpWT(RETA3,  WTZ,  EM) 

965  WTCHNG  =  WTX*WTY*WTZ 

PRIM  925,  BE  T  A3 ,  RF.TLOG,  SIGHA,  WTCHNG 
r  925  FORMAT ( Fl3 . 6 ,  Fit. 4,  Ft2.4.  F11.4) 

I F ( S  I GM A  -  0.0001  )  944 ,  930,  930 
*  930  CONTINUE 

Q  GO  To  944 

942  CALL  EXIT 
END 

0 


Q 


C 

c 

c 


0 

0 

0 

0 

0 

0 

o 


o 


FTN  1.4 
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c 
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SUBROUTINE  SIMPSON(SIMP) 

COMMON  SK(50)»  PROD(50),  A31,  A3?,  GET  A3 ,  XI3,  OMEGA,  EM 

NINCZ  s  40 

ZINC  -  0.025 

X ! 3  =  O.C 

CALL  8  I NOM ( 8  I N ■ : 

SIMP  =  Z  T  NC.-  (  3.  b  8 1  NT) 

BINTi  -  i.n/P-lNT 

PRINT  4.1,  R  E  T  A  3 ,  X  I  3  .  BINTI,  SIMP 

DO  30  JSIMP  =  1,  NINCZ 

AJSIMP  =  JSIMP 

X  I  3  s  ZINC* AJSIMP 

IFCJSIMP  -  NINCZ)  26,  27,  27 

26  IF<XMOnF{ JSIMP,  ?))  28,  29,  ?8 
2H  C At L  BIND M (BINT) 

panel  -  4,o#?inc/<bint*3.0) 
simp  =  SIMP  *  PANEL 
BINTi  =  1.0/BINT 

PRIM  4t,  B  E  T  A  .3  ,  X  I  3 ,  BINTI-  SIMP 
GO  TO  30 

29  CALL  PINOM(BINT) 

PANEL  =  ?.P*/IMC/(H[NT*3.0) 

SIMP  =  SIMP  +  PANEL 
R I N T 1  r  l.n/PINT 

PRIM  41,  13  E  T  A  3  ,  X  I  3 ,  B  I  N  T  t ,  SIMP 
GO  TO  30 

27  CALL  B I N O M ( R I N T  ) 

PANEL  s  7  I  M  C  /  (  P  I  N  T  *■  3 . 0  > 

SIMP  :  SIMP  ♦  PANEL 

BINTi  s  1 . P / A  I N T 

PRINT  41,  BETA  <, XI 3,  8  |  N I  1 ,  SIMP 

30  CONTINUE 

41  f  O  R  M  A  T  (  7  H  R  (:  T  A  3  s  ,  Flfl,*,  5*  NHXlJ  F7.4,  5X  lOMl.  /BlNT  a, 
1  El  2. 5.  MX  6MSIMP  3,  El  2. 5) 

RETURN 

END 


O 


c 


FTN  1,4  surroSti^Ix'?^!?.,  xr/j»  C212?7  ««  no.  t 

c  tTmuIJ°weora!"of  s°"kS  T#  BE  US6B  "•  S,KPS0N  INTtGRRTION 

IF  <  BE  T  A  -  0  .  0  0  0  3  >  U2,  \12,  1J3 

112  inxi  -  0,9)  114,  114,  115 

113  IF  <  BE  T  A  -  0.001  )  117#  117#  118 

117  I F  ( X I  -  0.8)  114,  H4.  H5 

118  IF(8ETA  -  0.003  )  122,  122,  123 

122  I F  <  X I  -  0.7)  114#  114,  115 


I F  <  BETA 
I F  <  X  I  - 
IF ( BETA 


0.01 


IFCXI  -  0.4) 

IF (BETA  -  0.03 


114#  114,  115 
)  127,  127,  128 
114,  114,  115 
)  132,  132,  115 
114,  114,  115 


inXI  -  0,0)  114,  114,  H5 

S  s  i.o 

RETURN 

CONTINUE 

PI  =  3.1415926536 
S  =  0.0 

DO  310  PM  =  1,  J00 

NPT  =  MM  -  1 

ANPT  a  NPT 

ANPH  =  ANPT  ♦  0.5 

ARG  =  RETA*ANPH*ANPH*PI*PI 

IF(ARG  -  175.0)  308,  3q8 ,  309 

FAC1  =  1.0/EXPFURG) 

GO  TO  306 
TERM  =  0.0 
GO  TO  331 

FAC2  =  COSF(  ANPfUPI  *x  I  ) 

I F ( XMODF ( NPT , 2 ) )  30*,  305,  304 

FAC3  =  2. 0/( ANPM*PI ) 

GO  TO  302 

F  AC3  =  -2. 0/( ANPMvP I ) 

TERM  s  FAC3  +  FAC1*MC2 

ir(ARSF(FAC2  )  -  1.0E-05'  303,  303,  327 
I F( NPT  -  3)  303#  303,  328 
IF(ARSF(TERM,$  ,  -  1.0F-05)  331,  331,  303 


IF ( A  R  S  F  ( TERM/ $  >  -  l.QF-05)  331  ,  331,  303 
S  =  S  +  TERM 
CONIT  IN'UE 
RETURN 

S  =  S  ♦  TERM 

FORNAT(6H  WT  =,  1.5.  MX.  SHS  ..  F«.S.  MX  6R-tRn  ,,  fel2.5) 
END 


r 


c. 


c 

C 

C 

/ 

c 

c 

c 

c 


c 


AT  021257 


PAGE  NO,  1 


PTN  1.4  DATE  08/25/67 

SUBROUTINE  PRQDSK 
C  SUBROUTINE  TO  DETERMINE  PRODUCT  OP  SK ( BET A1 ) *SK ( BET A2 ) 
COMMON  SK ( 50 ) »  PROD<50),  A31,  A32,  BET  A3,  XI3,  OMEGA,  EM 
I F (PROD ( 1 ) )  710.  711,  710 

710  RETURN 

711  BETA  =  BET  A3* A31*A31 

C  I F ( A 3 1  -  A32)  750  ,  751,  750  ' 

C  750  PRINT  735,  8ETA3,  BETA 

735  F ORM AT ( 9H  BETA  3  =,  Flo. 6,  10X,  8HBETA  1  r,  F10.6) 

C  GO  TO  782 
C  751  PRINT  733,  8FTA3,  BETA 

733  FORMAT ( 9H  BETA  3  =,  F10.6,  10X,  17HBETA  1  =  BETA  2  F10.6) 

C  782  PRINT  172 

172  F0KMAT(6H  XI  ,  5X ,  5HS  »  3X,  3HC  ,  40X,  5HSK(K>> 

CALL  SIMP1?(BETA) 


DO  720  MN  =  1,25 
7?0  PROD ( MN )  =  SK(MN) 

I F ( A 3 1  -  A32)  716,  717,  716 

716  BETA  =  BETA3*A32*A32 

C  PRINT  734,  0ETA3,  BETA 

734  FORMAT ( 9H  BETA  3  =,  FlO.6,  10X,  8HPETA  2  =,  F10.6) 
C  PRINT  172 

CALL  SIMP12(PETA) 

DO  730  MN  =  J .25 
730  PROD(MN)  =  PROD ( MN ) *SK ( MN ) 

RETURN 

717  DO  740  MN  =  1  ,25 

740  PROD(MN)  =  S« ( MN ) *SK { MN ) 

RETURN 


END 


\ 


FTN  1.4 


DATE  08/25/67  AT  021257 

SUBROUTINE  BINOM(BINT) 

COMMON  S  K ( 5  0 ) ,  BROD(50),  A31,  A32,  BET  A3 ,  X!3, 

BINT  =  1.0 

BCOEF  =  1.0 

DO  830  KB  I  =  1,  25 

K  :  KB  1 

A  KB  I  =  KBI 

BCOEF  =  BCOEF  *(FM-AKBI*1.0)/AKBI 
F AC  1  =  (OMEGA**(1  ,0/FM)  -  t.0)**K 
CALL  SxI(BFTA3,  XI 3,  S) 

F  AC2  s  S -  *  K 
CALI  PRCOSK 
F4CJ  =  PROD(KBI) 

TERM  =  BCOEF  ♦  F  A  C 1  *  F  A  C  2  *  F  AC3 
IF ( KR I  -  l )  832.  832,  828 
828  IM AHSF( TEBM/BINT)  -  l.oE-05)  831,  831,  832 
832  BINT  =  El  NT  *  TERM 
830  CONTINUE 


PAGE  NO. 
OMEGA,  EH 


1 


PE ( URN 

fl.vl  BINT  =  BINT  ♦  TERM 

840  FORM  A  T  (  4  H  K  13,  inx,  6H»lNT  =,  E12.5,  1 0  X  .  6HTERM  =,  H2.5) 
RETURN 
END 


c 


c 

i 

S 

* 

r 

w 

C 

c 

c 


FTN  1.4  DATE  08/25/67  AT  021257  PAGE  NO.  1 

SUBROUTINE  S I MP12 ( RET  A  ) 

C  SUBROUTINE  FOR  SIMPSON  INTEGRATION  OF  S(XU),  S(X12)  TO  GIVE  ALL 
COMMON  SK ( 50 ) »  PROD ( 50  ) ,  A31,  A32,  BETA3 ,  X  I  3 ,  OMEGA,  EM 
PI  =  3.1415926536 
XI  s  0.0 
ZINC  =  0.020 
NINCZ  =  50 

111  CALL  SXKBETA,  XI,  S> 

150  DO  152  MN  =  l,  25 
152  SK(MN)  =  Z I NC*S**MN/3 . 0 
C  s  1.0 

C  PRINT  171,  XI,  S,  C,  <SK(L>,  L  *  1.9) 

DO  180  JSIMP  =  1.  49 
AJSIMP  s  JSIMP 
XI  =  A JS I MP*7 I NC 
115  CALL  SXKBFTA,  XI,  S> 

139  IFtXMODFC JSIMP, 2) )  140,  141,  140 

140  DO  142  MN  =  I,  25 

PANEL  =  4,0*ZINOS**MN/3.0 
142  SKtMNl  =  SK(MN)  ♦  PANEL 
C  5  4.0 

C  PRINT  171,  XI,  S,  C,  (SKID,  L  =  1,9) 

GO  TO  180 

141  DO  144  MN  s  1,  25 

PANEL  =  2.0*ZlNC*S**MN/3.0 
144  SK(MN)  =  SK(MN)  ♦  PANEL 
C  =  2.0 

C  PRINT  171,  XI,  S,  C,  { SK  ( L )  »  L  =  1,9). 

171  FORM  A  T ( F  6 . 3 ,  F10.5,  F6.1,  9F10.5) 

180  CONTINUE 

C  PRINT  172,  ( S  K ( L i ,  L  =  1,  9) 

172  FORMAT ( 22H  1.000  0.00000  1.0,  9F10.5) 

C  PRINT  173 

173  FORMAT ( 1 H 1 ) 

RETURN 

END 


FIN  1,4  DATE  OR/25/67  AT  021257  PAGE  NO.  1 

SUBROUTINE  S I M  P  W  T ( B  E  T  A , W  T , F  M ) 

C  COMPUTATION  0  F  CONDUCTIVITY  CHANGE  AS  C  AI.CUL  A  I  H  D  «Y  WEIGHT 

C  CHANGE  FORMIJI  AS 

Pi  =  3 . 1 415926536 
sum  =  o.n 
DO  4f«0  K  =  1,  100 
AN  -  K  -  1 

TERM  ;  (1. n/(  AN  ♦0.5  )•*?)•  E  VPF(-0  E  TA*  P  I*PI*(  AN  ♦  0  .’>)**  2  ) 

IF(  TERM  -  1  .  "E'-OA)  481.  401  ,  4B0 
4  50  SUM  =  SUM  *  TW|,M 
4 H  l  W  T  s  ( ■?  .  0/ (  P  I  *•■"  I  )  )  *  (  SUM  ♦  T E MM) 

RETURN 

END 
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